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Periodic orbits in warped disks 
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Abstract. It is often assumed that a warped galaxy can be modeled by a set of rings. This paper verifies 
numerically the validity of this assumption by the study of periodic orbits populating a heavy self-gravitating 
warped disk. The phase space structure of a warped model reveals that the circular periodic orbits of a flat disk are 
transformed in quasi annular periodic orbits which conserve their stability. This lets us also explore the problem 
of the persistence of a large outer warp. In particular, the consistency of its orbits with the density distribution 
is checked as a function of the pattern speed. 
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1. Introduction 

In the conventional understanding of galaxies embedded 
in a kinematically hot spheroidal or triaxial halos of dark 
matter, tiie frequently observed pronounced outer warps 
of neutral hydrogen, interstellar dust and, to a lesser no- 
ticeable extent, optical disks set a challenge for explaining 
their presumably long-lived existence, and particularly, 
their frequent quasi-straight li ne of nodes (LON) , or pair 
of quasi-straight LON (Briggs 1990 ). At least half the spi- 
rals possess a detectable warped HI and stellar disk(e.g., 
Briggs 1990, and more should remain undetected due to 
projection confusion. 

Hypotheses like cosmic infall (Jang fc Bin ney |1999D , 



gravitational interactions (Hernquist [1991; Weinberg 



1995, 1998), normal modes (Sparke 1984; Sparke & Caser- 
tano 1988), misaligned dark halos (Kuijken 1991; Dubinski 
& Kuijken p^ ; Deba ttista fc Sellwood |l99|)' or magnetic 
fields (Battaner et al. 1990] ) have been proposed, without 
clearly satisfying all t he co nstraints pro vided by the obser- 
vations. See Binney ( |l992| ) or Kuijken ( ^000]) for reviews. 
But, as argued by Arnaboldi et al. ( 1997| ), Becquaert 



& Combes ( |1998[ ) and Reshetnikov & Combes ( |1998D , also 
in the context of polar rings, the evidence for non-self- 
gravitating HI disks is actually weak. Instead, the oppo- 
site assumption of coexistence of warps containing most of 
the dynamical mass is not in contradiction with the avail- 
able observations. Such an assumption of thick and flaring 
disks of visible and dark matter in the form of cold gas 
approximately proportional to neutral hydrogen has been 



proposed by Pfenniger et al. (1994) to account for many 
more known facts about spirals. 



For the warp problem, the self-gravitating disk as- 
sumption is attractive for several reasons. HI observations 
reveal that most of the galaxies possess an inner flat disk, 
while the war p dev elops beyond a specific radius (Briggs 
1990 ; Burton 1992). Such a comportment, which is the 



signature of a radical change of the dynamic along the ra- 
dius, can hardly be explained by a hot halo which would 
impose a uniform dynamic over a larger scale. In con- 
trast, a strong rotational support in the dark matter com- 
ponent offers a natural explanation for the straight lines 
of nodes. They would result as a natural consequence of 
the weak diffusive properties of angular momentum in a 
self-gravitating disk: if, foUowing for example an accretion 
event, angular momentum is deposited in a self-gravitating 
disk, bending its outer parts, the angular momentum will 
slowly diffuse across the whole disk in a monotonous way, 
because angular momentum is a quasi conserved quantity. 
As its diffusion is slow, the line of nodes must be quasi- 
straight and persist for many rotational periods. Finally, 
self-gravitating optical disks are in agreement with the of- 
ten stated maximum disk interpretation of observations of 
the Milky Way and other spirals (see for example Gerhard 

pooo| ). 
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In numerous models of warped galaxies, and also po- 
lar ring models, stars are assumed to move along circular 
rings which are tilted as a function of radius. Because stars 
are driven by periodic orbits, which are the backbone of 
galaxies, such models implicitly presuppose the existence 
of stable circular tilted periodic orbits. Yet, to our knowl- 
edge they have never been verified, and do not appear 
obvious in non-spherical geometries. It is far from obvious 
that a perturbation in the form of a warp will conserve 
the well known circular orbits of a flat disk. In order to 
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verify the latter assumptions, it is useful to find the exact 
periodic orbits existing in a warped disk. 

Moreover, understanding the stable periodic orbits al- 
lows us to grasp the basic properties of the other quasi- 
periodic orbits in an eflicient way, particularly when no 
analytical tools exist. Such an approach has been very 
useful in the context of barred galaxies for which the 
coniplexity of motion is substantial (e. g., Contopoulos & 
Papayannopoulos 198C| ; Pfenniger 1984). 

For studying the principal periodic orbits in a galactic 
potential, even rough representations of the potential are 
sufficient to yield the basic properties of the main peri- 
odic orbits, as that shown in numerous situations. In this 
work, the disk autogravitation will play a dominant role. 
Therefore, the foUowing results will be applicable as long 
as the local density is dominated by the disk. 

The alternative case would be the existence of a hot 
dark matter halo. In this case, one must admit a center 
dominated by the halo density, because the density of a 
warm gravitating system increases at the center. The fol- 
lowing implications should be explored: 

1) A hot halo is aligned with the central stellar disk. The 
outer gaseous disk acts as a tracer of the main periodic 
orbits in the halo. The only reason for the disk to warp 
is then to foUow the orbits generated by a vertical 1:1 
resonance. But in spheroidal halos, stable warped or- 
bits do not exist due to the invariance in azimuth of 
the potential. To avoid this problem, one has to take in 
account both the mass of the halo and the self-gravity 
of the disk (Sparke |1984| Sparke & Casertano |198SD . 
In triaxial halos more periodic 1:1 orbit families ex- 
ist. 1:1 resonances are strong o nly in rapidly rotating 
triaxial potentials (e. g. Mulder 1983 ), and associated 
orbits would be compatible with warps only outside 
the corotation radius: the dynamics woul d res emble 
the one of a rescaled triaxial bar (Pfenniger 1984). But 
such triaxial halos appear difhcult to reconcile with the 
axisymmetrization that foUows the inclusion of even a 
weak fraction of dissipative matter (Dubinski 1994 ). 

2) The gaseous disk is aligned with the outer part of a 
hot halo, and the optical disk is tilted. This case ap- 
pears dynamically highly unstable because the optical 
disks appear maximum (cf. Gerhard 2000), i. e., self- 
gravitating to a large extent. 

3) The halo foUows the warp of the disk. Such a halo 
can survive if it possesses enough angular momentum 
and therefore should be flattened. This case becomes 
very similar to the case of the massive disk, and the 
comportment of the periodic orbits is not changed dra- 
matically. 

In this paper we present succinctly a study of periodic 
orbits in a simple model of warped disk galaxies, in which 
the stability and consistency with the density distribution 
is also determined. The influence of a global rotation of 
the warp pattern is examined too. A subsequent study will 
investigate the orbits in self-consistent N-body models of 
warps. 



2. Model 

2.1. Warped disk model 

In order to compute periodic orbits, we have based our 
model on a superposition of three simple Miyamoto-Nagai 



potentials (Nagai & Miyamoto 1975). The three compo- 



nents can be viewed as representing respectively a bulge, 
a visible disk and a gas disk containing a large amount of 
dark matter. In cylindrical coordinates, the total potential 
can be written as : 



$o(i?,'^,z) = - 
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(a, + y 



(1) 



Observations suggest that warps seen on edge take the 
shape of an integral sign, i. e., the disk deformation is pro- 
portional to the cosine azimuthal angle (j) ^-nd increases 
with radius. Analytically, the warp is well represented by 
a term Az where : 



Az — wR cos (j), 



(2) 



and w is an adjustable parameter characterizing the 
warp's amplitude. 

Replacing z by z — Az in the potential (0) we get an 
analytical warped disk potential model : 



GM,, 



with 



di = y {z — wR^ cos < 



6? 



(3) 



(4) 



Note that the potential is bi-symmetric with respect to 
the y — Q plane, $(a;, y, z) — <&(— a;, y, —z). 

2.2. Parameters 

The parameters a^, bi and G Mi have been chosen to pro- 
vide a typical rotation curve, increasing linearly below 
3kpc and staying flat up to 40kpc, beyond which it de- 
creases smoothly. For the bulge, af, — O, bt = 1.5 and 
GMh = 0.066, which is equivalent to a Plummer sphere. 
The parameters for the visible and gas disks are respec- 
tively : ad = 6.5, ha = 0.5, GMd = 0.171, Ug = 25, 
bg = 0.5, GMg = 0.763. For convenience the length unit 
is chosen to be the kpc and time unit the Myr. Fixing the 
gravitational constant G to a value of 1, the mass unit cor- 
responds to 2.2 • 10^^Mq. With this set, the velocity unit 
is about 1000 km • s~^. Fig. |l] shows the isodensity curves 
corresponding to the model. Well above the galaxy disk, 
negative density regions do exist, however our results are 
not spoiled by negative densities as long as orbits do not 
cross such regions. The amplitude of the deformation w 
is chosen to be equal to 0.005 kpc~^, developing a warp 
of the same order as the one of the Ga laxy. For exam- 
ple, se e the values indicated by Burton ( 1992|) or Smart 
et al. ( 199^ ). Hence, the results obtained are applicable in 
typical cases of observed warped galaxies. 
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Fig. 1. 3D isodensity surface of the warped model at 5 ■ 
is sliced with respect to the y = O plane. 



10"^ and 5 • lO^^Mg • kpc"'' (w = 0.005 kpc^^). The surface 



2.3. Hamiltonian 

In the present model, we allow the possibility for a global 
rotation of the warped potential about the 2:-axis. This 
rotation is parametrised by the angular pattern speed Qp 
or, equivalently, by the corotating radius i?c, taking a neg- 
ative value for retrograde rotation. Thus, the Hamiltonian 
in the rotating frame of reference with angular speed Qp 
is : 

i? = i (pI+pI+pI) +^{x,y,z) - np{xpy-yp.^), (5) 

where the variablespa^, Py, andp^ are the respective canon- 
ical momenta of x, y, and z. The equations of motion to 
integrate are : 



y = Py - flpX 

z = Pz, 



with given by Eq. (||). 



Px 
Py 
Pz 



ripPy, 

^pPx ^ 



(6) 



2.4. N u merica I method 

The periodic orbits are found by numerically determining 
the fixed points of the 4D Poincare map (T) at y = O, 
a; < O generated by the equations of motion (see, e. g., 



where Ai and A2 represent a pair of reciprocal eigenvalues. 
A periodic orbit is stable only when bi and 62 are real and 
1 5i 1, 1 62 1 < 2. It is unstable in all other cases. If |&i| = 2 or 
1 62 1 = 2, or if |6i| = I62I, at least two eigenvalues are equal, 
Vr is degenerate and eventually allows a bifurcation. For 
a more complete description of the instability cases, see 
Pfenniger & Friedh ( |l993|) . 

In this work, we will focus our investigations on or- 
bits foUowing the potential disk. Hence, we choose initial 
conditions with a starting position on the line of nodes, 
here equivalent to the j/-axis. We concentratc on the main 
orbit families keeping the symmetry of the potential, thus 
the initial velocity is perpendicular to it (j)y{0) = 0) and 
points to arbitrary negative values of x (pi(0) < 0). In 
this way, the free parameters in the initial conditions are 
the Jacobi constant (w energy. H), the position y{0) along 
the y axis, and the velocity component Pz(0) along the z 
direction. This latter parameter will be small compared 
to PxiO), which avoids a situation where orbits stray too 
far from the disk and penetrate negative density regions. 
This point has been systematically checked because, a pri- 
ori orbits may explore regions far from the expected ones. 



3. Families of periodic orbits 



Pfenniger & Friedh 1993). The algorithm uses in particu- 3.1. Families witliout perturbation 



lar the method proposed by Henon (1982) giving the inter- 



section between an orbit and a surface, and a least squares 
stabilized Newton-Raphson root-finding procedure. 

The stability of the orbits are determined by the eigen- 
values of the Jacobian of the Poincare map (VT). Since 
the System is Hamiltonian and the motion is described by 
real numbers, the four eigenvalues of VT occur by conju- 
gate and inverse pairs and it is possible to condense the 
Information with two stability indexes 61 , 62 : 



-(A, + Ari), i = 1,2, 



(7) 



In order to understand the influence of the warp, we first 
look at the model without deformation {w = 0) which is 
completely axisymmetric. The rotation of the potential is 
set to zero (fip = 0). 

Below the orbit, families are shown with thcir initial 
starting point in the H—pz{0) diagram. In Fig.||, bottom, 
we recognize the circular orbit family (horizontal line at 
Pz{0) = 0). The stability indices for this family are traced 
at the top of Fig |[ In this particular case, bi indicates 
the stability in the galaxy plane while 62 corresponds to 
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Fig. 2. H-pz{0) phase space for u; = 0. 



the stability transverse to it. Both indices remain in the 
interval [—2, +2], insuring the stability for this family. 

At H = -0.0701, bi = +2, VT is also degenerate with 
two eigenvalues equal to —1. A bifurcation occurs through 
period doubling. This bifurcation coincides with the res- 
onance between the radial frequency k and the circular 
frequency fl (2k = 3 il). 

At H = -0.0748, -0.0566, -0.0451, -0.0345, and 
—0.0236, 62 — +2 generating transverse bifurcations 
through period doubling. These bifurcations coincide with 
the resonances between the transverse frequency i' and the 
circular frequency {2v = {2k + 1) fi, fc = 1, . . . 5). Orbit 
families from period doubling bifurcations will not be dis- 
cussed further in this paper. 

Transverse bifurcations keeping the same period oc- 
cur at = -0.0632, -0.0508, -0.0396, -0.0294, and 
—0.0158. In this case 62 = —2 and v — kO., k — 2, . . .Q. 
The two latter families have a non zero Pz(0) initial ve- 
locity and also oscillate round the z = O plane crossing it 
respectively 2k and 2{k + 1) times per period. AU these 
sub- families are marginally stable with 62 = —2 due to the 
axisymmetry of the potential. 



3.2. Influence ofthe warp 

The break in the z symmetry transforms the family of cir- 
cular orbits in the unperturbed model into a set of families 
pk (main families) with Pz(0) 7^ 0. The significance of the 
indice k will be explained further. Except at bifurcations, 
the pk families are stable. 

The shape of the corresponding orbits can be approx- 
imated with the foUowing parametric function: 

R{H) = i?o(i/)+^ii(i?)cos(20-7r), (8) 
z{H) = zo{H) + A,{H)cos<t>, (9) 

parametrised by the azimuthal angle cj). For increasing H, 
Ar increases from O to 1 kpc, while Az increases from O to 
8 kpc. This latter variation corresponds almost exactly to 
the amplitude of the warp as a function of radius [A^ = 
wR^). In other words, the main orbit family foUows the 
density maximum at any radius. 

Another influence of the warp is to move the bifur- 
cations toward slightly lower energies. The corresponding 
transverse bifurcations arise aX H = —0.0633, —0.0509, 
-0.0398, -0.0298, and -0.0176. Moreover, a new bifurca- 
tion takes place within the limits of the diagram at H ~ 
—0.0108, where the subfamily starting at iJ = —0.0176 
rcjoins the main family. The other families would also re- 
join the main family at higher energies, yet they would 
reach radii larger than y{0) = 50 kpc, which was chosen 
as a limit for our study of a warped disk. The subfamilies 
coming from the corresponding transverse bifurcations are 
given the symbol "sA:±" where k gives the frequency ratio 
v /n at the corresponding bifurcation and the sign is the 
one of the difference of pz(0) between the subfamily and 
the main family. These families evolve in the same way as 
family pk, i. e. they foUow the warped disk, so may be sig- 
nificantly populated in a real disk. Note that the families 
with an odd k conserve the bi-symmetry of the potential, 
which is not the case when k is even. Fig. ^ displays a set 
of three orbits extracted from different subfamilies, sA+, 
s5+ and s6— . The right panels show the elevation z of the 
orbits as a function of the azimuth (j) in comparison to the 
elevation ofthe corresponding main family (dashed lines). 
The more the energy grows, the more the subfamilies fol- 
low well the main orbit. Although in the absence of the 
warp all subfamilies are stable, this is no longer the case 
with a weak warp. This point will be discussed further 
below. 

3.3. Surfaces of section around bifurcations 

If the warp has only a slight influence on the shape of the 
orbits, it is important to look carefuUy at the behavior 
of the bifurcations in presence of the perturbation. Fig. |^ 
and Fig. || present a magnification of the H-pz{0) phase 
space around the first (rl) and second (r2) bifurcations. 
Here the periodic orbits are computed with w = 0.01 in 
order to increase the effect of the warp. The physical cor- 
respondent would be a galaxy with a warp twice as high 
as that observed in the Milky Way at a radius of 30 kpc. 
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Fig. 4. Projections over the y — O, z — O planes and the evolution of z as a function of azimuthal angle (f) of orbits 
of the families s4+ (H = —0.03), s5+ {H = —0.02) and s6— {H — —0.015). The dashed Une represents the potential 
minimum of the warped disk at a fixed radius which is similar to the main family. 



Except for the weak decreasing of Pz (0) with respect to 
H, the type of the first bifurcation (around H — —0.0633) 
remains similar to the one without perturbation (piteh- 
fork). The two new generated famihes (s2+ and s2— ) are 
both stable. They appear at the same energy as two other 
unstable famihes, namely z2+ and z2— . These are not 
seen in Fig. ^ because their projection merges with the 
p2 family. They are distinguished from it by a non zero 
z{0) and Py{Q)- Projections of families s2+ and s2— are 
plotted in Fig. §. 

In order to estimate the importance of orbits asso- 
ciated with a stable periodic orbit, we have computed 
the surface of the section {z, p z) near the bifurcation 
{H = —0.06). Moreover, it aUows us to know the num- 
ber of effective integrals of motion in this specific region 
of phase space. In Fig. 6 the stable p2 family is represented 
by the cross near the center (z = O, Pz = —0.0058) and 
is surrounded by quasi periodic families. The thinness of 
these curves indicates that the motion is well decoupled 



in z and a third integral exists. At Pz = 0.031, z = O 
and Pz — —0.043, z = O we find respectively the s2+ and 
s2— families. The upper and lower panels show details of 
the section around these two families. Very thin stability 
islands exist. For comparison, the width in pz of the sepa- 
ratrix (squares points) is about 1/100 times smaller than 
the width between families s2+ and s2— . This gives an in- 
dication of the phase space occupied by orbits associated 
with these families. The unstable z2+ and z2— families 
appear on the separatrix, at Pz = —0.0054, z — 0.427 and 
Pz = —0.0054, z = —0.427. Outside the separatrix, the 
invariant curves are regular until Pz = 0.06 where the mo- 
tion in z begins to be strongly coupled to the motion in x 
and y. 

The behavior of the second bifurcation r2 around 
H — —0.0509 is quite different. Fig. ^ reveals that the 
family s3+ is in fact the extension of the main family 
p2. This arises because the perturbation favors a higher 
Pz- The s3+ family is also stable. At H ^ -0.0508 the 
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Fig. 10. Shape of three among the four orbits hnked to the r2 bifurcation {H = —0.0509). The dot indicates the 
starting point of the computation and the arrow shows the direction of the corresponding initial velocity. 



families p3 and s3— are created simultaneously. They are 
respectively stable and unstable. The whole bifurcation 
forms a pitchfork with a symmetry breaking. The shapcs 
of three of the four families involvcd in this bifurcation are 
presented in Fig. The section aftcr the bifurcations is 
shown in Fig. I (i/ = -0.05). The stable p3 family appears 
at pz = —0.009 and is surrounded by a set of quasi peri- 
odic orbits embedded in the separatrix. The unstable s3— 
family belongs to the separatrix at Pz = —0.023. The s3+ 
family is located at Pz = 0.006 and its associated orbits 
occupy the crescent-shaped region between the separatrix. 
Outside, the invariant curves are regular until pz = —0.03 
where the motion in z is not longer decoupled. Looking at 
Fig. 6 and pl, one can expect a chaotic region around the 



separatrix. Due to the thinness of the region around the 
unstable points, it is very difficult to find it numerically. 

Our study of the bifurcation is limited to the two first 
bifurcations, rl and r2. The behavior of the foUowing bi- 
furcations are similar to the first two. The rk bifurcation 
with an odd (resp. even) k is of the same type as rl (resp. 
r2). This difference is related to the symmetry of the fam- 
ilies which depends on the parity of k. 

In summary, the warp does not destroy the strength of 
the main family, which at any H possesses a comfortable 
surrounding stable region. For an odd k, the effect of the 
warp is : (i) to change the stability of sk— subfamilies and 
(ii) to increase the phase space allowed to the associated 
sk+ families. For an even fc, the subfamilies are quasi un- 
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changed, preserving a very thin phase space around the 
stable subfamiHes skzL. 



3.4. Influence of rotation 

To study the influence of the rotation, we first look for 
Lagrangian points. The latter are found by setting the 
right terms of equations to zero. Substituting p^, Py 
and pz, the Lagrangian points are solutions of the foUow- 
ing equations : 



(10) 
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Fig. 6. Section {z,pz) a± H — —0.06, near the bifurcation 
rl. The crosses represent the three stable periodic families 
s2+ (top), p2 (middle) and s2— (bottom) and the two 
unstable szA- (right) and sz— (left). The upper and lower 
panels show a zoom of the section around s2+ and s2— . 
In this panels the separatrix is marked with squares while 
in the middle it is marked with a bold line. The crosses 
indicate the position of the periodic families. 

a.$.„ = 0. 

With the form (||) of the warped potential, one can easily 
see that the z coordinate must satisfy : 



z = Az, 

with Az defined in equation 
equations of (|o[) to become 



(11) 

This causes the two first 
(12) 



which corresponds to the dcfinition of the corotation, 
where $0 is the potential without pcrturbation (Eq. (0)). 
Thus the Lagrangian points are degenerated, forming an 
annulus foUowing the density maximum with a projection 
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on the z = O plane corresponding to the corotation with- 
out perturbation (w = 0). 

Despite the fact that a direct and retrograde rotations 
respectively add and remove resonances, the type of bi- 
furcations as well as the stabihty and the shape of orbits 
are not affected by a global pattern speed. This has been 
tested in the range —30 > Rc and Rc > 30 kpc, where Rc 
is the corotation radius and a negative value corresponds 
to a retrograde rotation. 

For a radius less than 40kpc, both prograde and ret- 
rograde 1:1 resonances are not observable in this range of 
pattern speed. For a corotation of 7 kpc they appear only 
beyond 40 kpc (resp. 30 kpc). Thus they have not been 
studied. 

The main influence appears when we look at the consis- 
tency of orbits with the mass density. If the warp is mostly 
self-gravitating and made of thin and distinct tube orbits, 
one can check the self-consistency constraint by noting 



that the spatial occupation of a periodic orbit is locally 
inversely proportional to its local speed {tso ^ l/l'^D- The 
reason is explained by the fact that at a given point of the 
orbit, the local speed remains constant in time. This ar- 
gument is not available at exceptional points, for example 
the points where an orbit crosses itself. Since the density 
is proportional to the spatial occupation time, it must 
also be proportional to the inverse speed along the orbit 
{p ~ tso ~ l/l"*^!)- Strictly this condition is only fulfilled 
by a structure entirely made of distinct exactly periodic 
orbits, such as a disk made of circular orbits. Nevertheless 
the check is useful in this problem because few hot orbits 
far from periodic round orbits are expected to exist. 

In practice, the consistency has been calculated using 
the indice / defined by the norm : 



1 



/(f7„iJ) = — 



2n 



2tt 



1/2 



with 



Ap(0) = 



Au((/)) 



p{(f)) - pmin 



Pn 



(13) 



(14) 



(15) 



and u(0) = l/\v{4>)\- The values pmin, Pmax, u^in and 
Umax are calculated over a period of the azimuthal angle 
0. When Ap and Au vary in the same direction, / = O 
and the orbit is consistent. If they vary inversely, as two 
cosines with a phase difference of tt, / = ^/^/271. 

We have tested the consistency of families pk with re- 
spect to the energy (Fig. ^l]) and radius (Fig. ^2|), which is 
more convenient in galactic dynamics. This as been calcu- 
lated for different (direct and retrograde) pattern speeds 
between —30 > Rc and Rc > 30 kpc corresponding to 
-4.5-10"^ <ilp < 4.5-10"^ Myr^\ The white regions cor- 
respond to / = O (consistency) while the darker gray cor- 
respond to a value of ^JH jl-K (inconsistency) . The shaded 
parts give the limit of the computation, either because of 
the corotation or because of the positive energy regions. 
The black pattern points out the missing data due to com- 
putational difficulty arising because of the proximity of the 
forbidden regions. 

For a direct rotation, except in small regions near the 
corotation, the inverse local speed along an orbit varyies 
exactly in opposition to the density (/ = \/7r/27r). This 
causes it to depopulate the higher density regions to the 
advantage of the lower. The density distribution is also 
slowly modified. However, for a slowly retrograde rotation 
a zone appears around 3 kpc, where the inverse local speed 
varies exactly as the density. This latter is also reinforced 
and the potential is self-consistent. This zone grows with 
increasing rotation and a second zone appears for larger 
radii. For a corotation smaller than Rc = —50 kpc the 
zones of consistency cover the whole disk under the curve 
of zero energy. 
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Fig. 11. Consistency of orbit occupation with the mass 
density as a function of the energy and the global rota- 
tion. The dark gray corresponds to inconsistent regions 
while the white are consistent regions. The black parts 
are missing data. The shading part represents the limit of 
the corotation. 



a o 




Fig. 12. Same graph than in Fig. |Tl| but as a function of 
radius R. The upper shading part represents the limit of 
the corotation. While the lower is the region where H > 0. 
The vertical resonances {ly — k fl, k — 2, . . . 10) are drawn 
in solid lines. 



The same work applied to the subfamilies sk± reveals 
that they are clearly inconsistent with the density distri- 
bution, but these familics are also less rclcvant with the 
assumption of almost circular rotation. 

4. Conciusions 

The previous description of periodic orbits in a warped 
analytical potential allows us to draw the following con- 
ciusions. 



1) The circular orbit family with an axisymmetric poten- 
tial survives the warp and becomes quasi-annular with 
a tilt corresponding to the warped galactic disk. Thus, 
this work confirms the intuitive idea that the orbits 
foUow the warp. Without this confirmation, the ring 
model of warps would have remained baseless. Note 
that the ring model has been used in the past in dif- 
ferent contexts, but if the underlying circular orbits do 
not exist, the ring model is nonphysical. 

2) The existing transverse resonances in the axisymmetric 
case are also preserved in the warped case but occur 
at slightly lower energies. Subfamilies with the same 
period as the quasi-annular family are generated and 
most of them keep their stability. In our model, the 
perturbation induced by the warp is not strong enough 
to generate observable chaotic regions. 

3) Since no bifurcation with v = ± Qp is generated by 
the warp, it is not possible in our model to excite a 
1:1 normal mode, which is often invoked to explain 
the galactic warps. Other modes are not favored either 
because the stable regions around their corresponding 
subfamilies remain very thin. 

4) The study of the consistency of the density distribution 
underlines the importance of the global rotation. A di- 
rect rotation of the warp pattern generates orbits de- 
stroying the density distribution, but on the contrary, 
orbits can locally reinforce the density in the presence 
of a retrograde rotation. This could help maintaining 
external warped regions. 

The behavior of the discussed periodic orbits will be 
examined in N-body models and the results will be the 
object of a future paper. 
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